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Thermally-densified hafnium terephthalate UiO-66(Hf) is 
shown to exhibit the strongest isotropic negative thermal 
expansion (NTE) effect yet reported for a metal-organic 
framework (MOF). Incorporation of correlated vacancy 
defects within the framework affects both the extent of 
thermal densification and the magnitude of NTE observed 
in the densified product. We thus demonstrate that defect 
inclusion can be used to tune systematically the physical 
behaviour of a MOF. 

Introduction 

It has long been recognised that defects and their correlations 
play a central role in determining the physical properties of 
many functional materials: the mass transport pathways of 
fast-ion conductors, 1 2 charge localisation in high-temperature 
superconductors, 3 4 and emergence of polar nanoregions in re- 
laxor ferroelectrics 5 are just three examples. The recent dis¬ 
covery that correlated defects can be systematically introduced 
into the structures of some canonical metal-organic frame¬ 
works (MOFs) 6 suggests the possibility of establishing de¬ 
fect/property relationships in this broad family of materials 
long-favoured for its chemical versatility. 7 9 Perhaps the most 
obvious benefit of determining such relationships would be 
the ability to exploit defects and non-stoichiometry to develop 
MOFs with particularly attractive or otherwise-inaccessible 
physical properties. 

The few recent studies of defective MOFs have focussed ei¬ 
ther on the nature of defect inclusion itself or on how defects 
influence chemical properties such as catalytic activity and 
guest sorption. 10-13 p[ ere 

we ask a different question: namely, 


f Electronic Supplementary Information (ESI) available: Synthesis, experi¬ 
mental methods, and sample characterisation; X-ray powder diffraction re¬ 
finement details. 
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is it possible that defects can influence the physical proper¬ 
ties of MOFs? In attempting to answer this question, we 
focus in this proof-of-principle study on the thermomechan¬ 
ical behaviour of the canonical hafnium terephthalate MOF 
UiO-66(Hf). 14 Our choice of this specific pairing of physical 
property and chemical system is motivated by the following 
considerations. First, the response of framework geometry to 
changes in temperature is at once both straightforwardly mea¬ 
sured and strongly characteristic of the underlying elastic be¬ 
haviour of the material in question. 15 16 Second, the empirical 
propensity for MOFs to exhibit anomalous mechanical effects 
such as thermal contraction (i.e. negative thermal expansion, 
NTE; Refs. fT7ifl9l ) has led to a general expectation that such 
properties will be especially sensitive to the existence and na¬ 
ture of structural defects. 9 And third, UiO-66(Hf) is the MOF 
for which there is arguably the greatest and best-understood 
chemical control over defect incorporation J 6 ! 1 

The idealised structure of non-defective UiO-66(Hf) 
consists of a face-centred cubic array of oxyhydroxy- 
hafnium(IV) clusters connected via terephthalate (benzene- 
1,4-dicarboxylate, bdc 2- ) linkers [Fig. [lJa)]J® The particular 
defects we study here are introduced by substituting formate 
for terephthalate during framework synthesis. As a mono¬ 
topic (capping) ligand, formate reduces the connectivity of 
the framework. This, in turn, encourages vacancies of en¬ 
tire hafnium clusters. Formate content can be increased at 
least as far as the situation where, of the 12 carboxylate lig¬ 
ands of each hafnium cluster, four ligands are formates and 
eight are terephthalate carboxylates [Fig. |TJb)]. Ed While this 
substitution results in substantive structural modification over 
the unit-cell (nanometre) length scale—indeed one in four 
hafnium clusters are vacant in this case—the actual coordi¬ 
nation environment within the clusters is essentially indistin¬ 
guishable in the defective and defect-free modifications.® 22 23 
In order to amplify the structural consequences of defect in¬ 
corporation we exploit the high-temperature ligand elimina¬ 
tion reaction known to occur in this system. 11 24 25 Whereas 
this step involves elimination of H 2 O for defect-free UiO-66, 
the lability of non-bridging formate ligands is such that it is 
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pendent on defect concentration. Quantum mechanical cal¬ 
culations help interpret the microscopic origin of this volume 
reduction. Second, all densified-UiO-66(Hf) samples exhibit 
isotropic NTE behaviour that is many times stronger than that 
of any other MOF, and in one case even satisfies the condi¬ 
tion of “colossal” NTE (its volumetric coefficient of thermal 
expansion ay = dV/VdT ~ —100MK -1 ). This is the first 
report of colossal NTE in an isotropic MOF. Third, the mag¬ 
nitude of NTE can also be tuned systematically by defect con¬ 
centration. Our paper concludes with a discussion of the likely 
mechanisms by which defect inclusion might influence elastic 
properties of this particular MOF, the challenges posed for un¬ 
derstanding NTE in defective materials, and the implications 
for exploiting defect/property relationships in other MOF sys¬ 
tems. 

Results and discussion 



Fig. 1 The ideal non-defective UiO-66(Hf) consists of 
oxyhydroxyhafnium(IV) clusters connected via terephthalate linkers 
into a face-centred cubic feu topology network. A unit cell is shown 
in panel (a), with the Hf clusters represented as blue truncated 
octahedra and linkers as red rectangles. The atomic-scale structure 
of a single cluster is shown in detail in the top-left of the panel, (b) 
at high concentrations of formic acid promotes the formation of 
nanodomains of the primitive reo topology, which contains 
correlated linker and hafnium cluster absences. Absent hafnium 
clusters are shown in skeletal form. In spite of the differences in 
framework topology, the cluster coordination environments of feu 
(c) and reo (d) UiO-66(Hf) are very similar.®! After ligand 
elimination (right-hand panels of (c),(d)) the two cluster types are 
distinct in both coordination and symmetry. Atoms are coloured as 
follows: Hf, blue; O, red; and C, grey. Atoms lost during thermal 
elimination are coloured black. 


formic acid that is eliminated in our defective UiO-66, result¬ 
ing in markedly different cluster geometries for the two cases 
[Fig.0c,d)].™ 

Making use of in situ variable-temperature synchrotron X- 
ray diffraction measurements of a series of UiO-66(Hf) sam¬ 
ples with a range of defect concentrations, we proceed to 
demonstrate the following. First, ligand elimination results 
in framework densification, the magnitude of which is de¬ 


Having used the method of Ref. [6 to prepare a series of six 
UiO-66(Hf) samples with a range of defect concentrations, we 
sought first to establish that our samples included defects of 
the type discussed above. X-ray powder diffraction measure¬ 
ments revealed the existence of structured diffuse scattering in 
the form of broad superlattice reflections that are symmetry- 
forbidden in the defect-free system but are associated with 
nanodomains of correlated cluster vacancies in defective UiO- 
66 [Fig. [2ja)]. The ratio of the intensity of this diffuse scat¬ 
tering to that of the parent reflections is a measurable quantity 
that is sensitive to defect concentration. We find an approxi¬ 
mately linear relationship between this ratio and the quantity 
of terephthalic acid used during synthesis, suggesting that our 
ensemble of defective UiO-66(Hf) samples spans a suitable 
range of defect concentrations [Fig.|2fa,b)]. 

The onset temperature of ligand elimination was deter¬ 
mined using thermogravimetric analysis (TGA). The TGA 
traces of all our samples are qualitatively similar and can 
be i nterpret ed on the basis of previous studies of this sys¬ 
tem: 112527 each reveals three primary stages of mass loss, 
corresponding in turn to adsorbate volatilisation (^100 °C), 
ligand elimination (250-350 °C) and framework decomposi¬ 
tion (550-600°C) [Figs.gc), ESI Figs. l&2f]. Qualitative 
trends in defect concentrations can be inferred from the rel¬ 
ative mass losses associated with the second and third stages 
[Am e ii m and Am^ ecomp , respectively, in Fig. |2|c,d)]. Frame¬ 
work decomposition corresponds to the loss of terephthalate 
(leaving behind Hf 02 ) and so larger values of Amdecomp are 
expected for the least defective examples, which should con¬ 
tain the highest relative concentrations of terephthalate. This 
is precisely what we observe [Fig.[2jd)]. Fikewise the value of 
Am e iim should decrease as the defect concentration is reduced 
(since H 2 O is eliminated instead of the heavier HCOOH), and 
again this trend is observed in practice. So both X-ray diffrac- 
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Fig. 2 (a) A diffraction pattern of the most defective sample of 
UiO-66(Hf). Structure diffuse scattering located at primitive 
superlattice positions is a signature of the presence of correlated 
defects in UiO-66(Hf). The most intense of these broad reflections 
are the (100) and (110). On increasing the quantity of H 2 bdc in the 
reaction mixture, the concentrations of defects decreases, and so the 
relative intensity of these primitive reflections decreases ((a) inset 
and (b)). (c) The TGA trace of UiO-66(Hf) shows three steps 
corresponding to, in turn, guest volatilisation, ligand elimination and 
framework decomposition, (d) The decrease in defect concentration 
is evident in the variation with bdc concentration of normalised 
mass losses both on ligand elimination which decreases, 

and framework decomposition (Aradecamp)- Fitted curves are 
included as guides to the eye. 


loss (i.e. up to 250 °C) the most noticeable changes to the 
diffraction pattern were in the relative intensities of reflections 
sensitive to solvent occupation within the framework pores. 
There was only a very small change (~0.1 %) in the size of the 
unit cell within this regime. These observations are consis¬ 
tent with our inte rpreta tion of the TGA trace and with previ¬ 
ous investigations l n 25 In contrast, the ligand elimination pro¬ 
cess that occurs between 250 and 350 °C resulted in an order- 
of-magnitude larger decrease in the cubic unit cell parameter 
that has not previously been observed in defect-free samples 
[Fig.|3ja)]. The diffraction pattern of the densified framework 
formed at this higher temperature can still be accounted for by 
a model based on the ambient-temperature structure, adjusted 
for the change in lattice dimensions; however the scattering in¬ 
tensity in high-angle Bragg peaks was much reduced such that 
detailed structural refinement did not prove possible. On sub¬ 
sequent cooling from 340 to 100 °C all samples showed NTE 
[Fig. [3ja,b)] . The two processes of densification and NTE, 
together with their dependence on defect concentration, are 
discussed in turn below. 


Thermal densification 

In the absence of a diffraction-based structural model for the 
densified framework, we turned to quantum mechanical cal¬ 
culations to verify that the structural changes associated with 
formate elimination were consistent with the magnitude of 
volume collapse observed experimentally. From a computa¬ 
tional viewpoint it would have been unfeasible to use ab ini¬ 
tio methods to relax large atomistic configurations describing 


tion and TGA measurements point to a systematic variation in 
defect concentration across our series of UiO-66(Hf) samples. 

In order to determine the effect of defect concentra¬ 
tion on thermomechanical behaviour, we measured variable- 
temperature X-ray powder diffraction patterns for all six sam¬ 
ples. Because we were interested in determining the thermal 
expansivity of the phases produced by ligand elimination, we 
paid particular care to ensuring that our samples were heated 
sufficiently to drive this elimination step but not so far as to 
initiate framework decomposition. In practice this dictated an 
experimental protocol of heating each sample to a temperature 
of 340 °C and then measuring thermal expansivity on subse¬ 
quent cooling to 100 °C. So the temperature range probed in 
these experiments corresponds to the first two stages observed 
in the TGA measurements described above. 

For all samples the same qualitative behaviour was ob¬ 
served (although, as will form the focus of subsequent dis¬ 
cussion, the magnitudes of these changes varied from sam¬ 
ple to sample). On heating through the first stage of mass 



T (°C) 2 0(°) 


Fig. 3 (a) Evolution of the lattice parameter of the UiO-66(Hf) 
sample synthesised with the largest quantity of F^bdc (5.0 mmol) - 
i.e. the lowest defect concentration — on heating from room 
temperature to 340 °C then cooling to 100 °C. A linear fit to the 
lattice parameter over the temperature range 340-160 °C is shown, 
for the determination of ay. (b) The evolution of the most intense 
peak, the (111) reflection, throughout this process. 
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the nanodomain defect structure reported in Ref. 6j Instead 
we focussed on an ordered defect structure corresponding to 
substitution of four equatorial carboxylate sites by formate on 
each hafnium cluster; this is the ‘ordered-reo’ model of Ref. 
|6 For this model, ligand elimination corresponds to the chem¬ 
ical process 

Hf 6 0 4 (0H)4(bdc)4(00CH)4 Hf 6 0 8 (bdc) 4 + 4HCOOHt. 

In silico relaxation of structural models for the framework 
compositions represented on either side of this reaction equa¬ 
tion suggested that a 1.0% change in lattice parameters 
is expected upon HCOOH elimination: a = 20.862 A for 
Hf 6 0 4 (0H)4(bdc) 4 (00CH)4 and 20.652 A for Hf 6 O 8 (bdc) 4 0 
In contrast, it is already known that elimination of water from 
defect-free UiO-66 does not result in any significant volume 
change. 26 By inspection of the relaxed structural model for 
Hf 60 8 (bdc) 4 , we conclude that distortions of the Hf6 octahe- 
dra are responsible for this framework densification. Elimina¬ 
tion of HCOOH amplifies the tetragonal distortion of Hf6 oc- 
tahedra, reducing the effective octahedral volume by 4.2% and 
bending the cluster—bdc—cluster connection so as to couple 
local distortion to macroscopic volume contraction [ESI Fig. 4 
f]. The relaxed structural models are available as crystallo¬ 
graphic information files as part of the ESI.f 

If, as this analysis would suggest, the process of frame¬ 
work densification depends on the presence of structural de¬ 
fects, then we should expect stronger densification in UiO- 
66(Hf) samples with larger defect concentrations. By heat¬ 
ing our samples to sufficiently-high temperatures to ensure 
complete densification (and working within the constraints 
of available synchrotron beam time), we found this ex¬ 
pected property/defect relationship to be borne out in practice 
[Fig. |4|a,b)]. Our most defective samples showed densifica- 
tions that were larger by ~1% (of the total lattice parameter) 
than for the least defective samples; this result is consistent 
not only with the results of our ab initio calculations but also 
with previous experimental reports. 26 

Negative thermal expansion 

Perhaps the most intriguing physical property of densified- 
UiO-66(Hf) is its NTE response, which is evident in the in¬ 
crease of the cubic unit cell parameter on cooling from 340 °C 
[Fig.[3ja)] and in the anomalous thermal shifts of Bragg peaks 
in the raw diffraction data themselves [Fig.^b)]. While NTE 
is increasingly frequently observed amongst MOFs, there are 
two particular aspects of the behaviour we observe here that 


t The tendency for calculations to overestimate lattice parameters for organic 
and organic-containing materials is well-known phenomenon and arises due 
to the difficulty of correctly accounting for dispersive interactions. It does not 
impact the analysis of structural and lattice parameter trends in families of 
materials of similar composition and porosity . 28 
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Fig. 4 (a) The evolution of the lattice parameter through the ligand 
elimination process for a series of samples with different defect 
concentrations. To correct for minor variations in the onset 
temperature of ligand elimination between sample, the curves are 
offset in temperature such that the maximum change in volume 
corresponds to 300 °C. (b) The densification (percentage linear 
decrease) after ligand elimination, measured at 380 °C, highlighted 
in (a) by a dashed grey line, (c) Coefficients of thermal expansion 
for isotropic NTE MOFs. The two values given for 
densified-UiO-66(Hf) correspond to upper and lower defect 
concentrations.h® 9 (d) Variation in NTE behaviour as a function of 
decreasing defect concentration (increasing quantity of H 2 bdc in the 
reaction mixture). 


are especially attractive: on the one hand, the cubic crys¬ 
tal symmetry of defective UiO-66 means that its NTE is 
isotropic— i.e., it occurs with constant magnitude in all crys¬ 
tal directions; and, on the other hand, the magnitude itself is 
really very large. Thermal expansion effects in different ma¬ 
terials can be compared via the lattice expansivities 

with typical (positive) thermal expansion of engineering ma¬ 
terials corresponding to values of ay in the range 10- 
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30MK -1 . 29 Here, the lattice parameter data of Fig.^a) cor¬ 
respond to a thermal expansivity of ay = —97(4) MK -1 over 
the temperature range 160-340 °C, which is comparable in 
magnitude to that of “colossal” thermal expansion materials 
(|a| ~ 100MK” 1 ) such as Ag 3 [Co(CN) 6 ]PAlthough we ob¬ 
serve an even more rapid increase in lattice parameter on fur¬ 
ther cooling from 160-100 °C, we cannot rule out the pos¬ 
sibility that this feature may arise from re-uptake of guest 
molecules at these lower temperatures, and hence do not at¬ 
tach particular weight to that observation here. We have also 
carried out separate variable temperature diffraction measure¬ 
ments confirming the reversibility of the NTE behaviour [ESI 
Fig. 4]. 

Isotropic NTE behaviour of this magnitude is rare indeed. 
First discovered in ZrW20s in the late 1960s, 31 isotropic NTE 
is now known to occur in a handful of materials, such as 
ZrW 2 C >8 itself and related substitutional variants, 31 33 the per- 
ovskite analogues ScF 3 and Re0 3 BUS the cuprite-structured 
oxides Ag 20 and CU 2 O , 37 the Zn x Cdi_ x (CN )2 family and an¬ 
hydrous Prussian Blue analogues, 38 41 and the three MOFs 
HKUST-1, MOF-5, and MOF-14. EH Recent theoretical cal¬ 
culations suggest that isotropic NTE may occur in a wider 
range of MOF materials that currently experimentally reported 
and re-emphasise the important role of topology. 42 Of these 
various systems, single-network Cd(CN )2 is the only colossal 
NTE material (ay = — 100.5(15) MK _1 )J 41 Despite decades 
of theoretical and experimental studies there is only general 
consensus that low-energy transverse vibrational modes play 
at least some role in the NTE mechanisms for most of these 
materials. The details of NTE mechanisms remain contro¬ 
versial, with new studies frequently demanding reinterpreta¬ 
tion of long-held assumptions as and when they appear. 17 43 45 
Here we avoid speculating on a detailed NTE mechanism and 
focus instead on the measurement and defect dependence of 
the phenomenon. We note simply that, with its open frame¬ 
work structure of high-nuclearity clusters connected by long, 
light, and flexible organic linkers, densified-UiO-66(Hf) cer¬ 
tainly shares many of the basic design features thought to 
favour low-energy NTE modes. Fig.|4|d) places the isotropic 
NTE behaviour of UiO-66(Hf) in the context of the three other 
MOFs known to show the effect. 17 19 

As anomalous as the thermal expansion behaviour of 
densified-UiO-66(Hf) is, our key interest is in the possibility 
that the magnitude of NTE might be sensitive to defect con¬ 
centration. What sensitivity, if any, do we actually expect? An 
earlier study of the geometric contribution to NTE in molec¬ 
ular frameworks suggested that defects of the type thought to 
occur in UiO-66(Hf) should actually encourage NTE because 
they act to increase the relative contribution of NTE vibra¬ 
tional modes to the overall framework dynamics.®! Likewise, 
even simple Maxwellian counting suggests that the balance 
of degrees of freedom and geometric constraints would be 


tipped increasingly in favour of the former (and hence NTE) 
as linkers are removed from a three-dimensional network. 47 
Not taken into account in either analysis, however, is the ob¬ 
servation we make here that defects are also linked to frame¬ 
work densification. Empirically, one finds that higher den¬ 
sities usually result in weaker NTE—this is exactly what is 
observed in zeolites, Prussian blues, and the Zn x Cdi_ x (CN )2 
family alike. 38 40 48 To complicate further this balance of com¬ 
peting effects, we remark that the ligand elimination process 
alters the point symmetry of the Hf6 clusters in a way that 
may itself affect the dynamics—and hence the NTE—of de¬ 
fective UiO-66(Hf). The point symmetries of both defect-free 
and defective Hf6 clusters (D^d and D^, respectively) are in¬ 
compatible with one another, as they are with the point sym¬ 
metry of the crystallographic site on which they both sit (Oh) 
[Fig.Ec ,d)]. Consequently, the effective energy potential gov¬ 
erning transverse vibrational motion—which, by analogy to 
better-understood NTE materials, is likely to be implicated in 
NTE—will not be constrained by symmetry and hence will 
vary from site to site throughout the defective lattice. De¬ 
spite the differences in local symmetry between the reo nan¬ 
odomains and feu matrix, both phases retain their cubic sym¬ 
metry on the unit cell level. The thermal expansion within 
each domain will therefore remain isotropic. Differences in 
thermal expansion between phases will lead to the develop¬ 
ment of thermal strains at the domain walls, which would in 
principle manifest themselves as differential peak broadening 
between primitive and face-centered reflections on heating. 
The large strains produced by the densification precluded our 
measurement of this comparatively subtle effect. The effect 
of defect concentration on the average effective potential, its 
asymmetry, and its variance, collectively poses an interesting 
challenge for computational and theoretical methods of under¬ 
standing NTE lattice dynamics in this system. 

What we observe in practice is that our UiO-66(Hf) samples 
with higher defect concentrations actually show reduced NTE 
effects relative to defect-poor samples. The coefficients of 
thermal expansion determined from our variable-temperature 
X-ray powder diffraction measurements range systematically 
from —70(2) to —97(4)MK _1 as defect concentrations are 
reduced [Fig.j4jc)]. So any effect of increased flexibility aris¬ 
ing from the reduced framework connectivity at high defect 
concentrations must be outweighed by considerations of den¬ 
sification and/or local symmetry. Whatever the microscopic 
origin for this defect dependence, we note that the chemical 
control over NTE behaviour we observe here would usually 
be associated with aliovalent substitution, as exploited in con¬ 
ventional solid-state NTE compounds (e.g. the ZrV x P 2- x 07 
family- 2 ). That we can achieve similar control over the phys¬ 
ical properties of MOFs through the use of defect chemistry 
is precisely the proof-of-principle needed to demonstrate that 
defects might be used as a design element in their own right 
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when engineering MOFs with targeted properties. 


Concluding remarks 

One of the key motivations for exploring NTE behaviour in de¬ 
fective MOFs was always that NTE is very often diagnostic of 
other lattice-dynamical anomalies, such as pressure-induced 
amorphisation and/or mechanical softening. Consequently the 
variation in magnitude of NTE documented in this study im¬ 
plies that changes in defect concentrations are likely also to 
affect the mechanical stability of MOFs in high pressure en¬ 
vironments and their resistance to amorphisation. We find 
that this relationship is likely to be entirely counterintuitive 
for defective UiOs (at the very least): reduced NTE nearly 
always translates to greater pressure stability, and so higher 
defect concentrations might now be expected actually to frus¬ 
trate amorphisation rather than to encourage it. In turn this 
may provide a mechanism to improve the mechanical proper¬ 
ties of MOFs, such as their machinability and durability un¬ 
der operating conditions. Whether this bizarre relationship of 
defect-driven toughening actually holds in practice remains to 
be shown, but what is clear even at this stage is that further 
studies of defect/property relationships in MOFs are set to 
challenge our collective intuition of the physical consequences 
of structural disorder. 

While we have focussed heavily on the implications of vari¬ 
able NTE in the UiO -66 system, we note for completeness that 
the control over framework densification evident in Fig. |4|b) 
may highlight another means of exploiting defects in MOFs. It 
has recently been shown that densification of guest-containing 
MOFs through pressure or mechanical grinding can be an ef¬ 
fective method of immobilising dangerous species ( e.g . ra¬ 
dioactive iodine ). 4950 Perhaps the incorporation of defects 
might simultaneously increase capacity for adsorption of these 
guest molecules while also trapping these molecules more ef¬ 
ficiently as a result of an increased density change on frame¬ 
work collapse. 

In summary, we have shown for the first time how the ther¬ 
momechanical properties of a MOF can be systematically var¬ 
ied via controlled incorporation of defects. The thermome¬ 
chanical properties we observe in defective UiO- 66 (Hf) are 
extremely unusual, and include the strongest isotropic NTE 
effect ever observed for a MOF. This property in itself may 
have application in counteracting the positive thermal expan¬ 
sion behaviour of other engineering materials. Yet the most 
important result remains a conceptual one: namely that de¬ 
fects might now realistically play a key role in the design of 
new classes of functional MOFs. 


Methods 

Preparation of defective UiO-66(Hf) 

We synthesised defective UiO- 66 (Hf) using the method de¬ 
scribed for large scale synthesis in Ref. 6 . The most defective 
sample was synthesised as follows. HfCU (3 mmol, 99.9% 
(metals basis, <0.5% Zr), Alfa Aesar) and terephthalic acid 
(2.5 mmol, Sigma Aldrich) were added to a 250 ml Schott bot¬ 
tle and dissolved in N,N dimethylformamide (DMF) (40 ml, 
Sigma Aldrich). This mixture was sonicated until the metal 
and ligand had dissolved and then formic acid (20 ml, 95%, 
Sigma Aldrich) was added. The bottle was sealed then placed 
in an oven at 120 °C for 24 h. The product was isolated by 
filtration, washed (DMF) and then heated at 60 °C with DMF 
for 3 days to ensure the removal any residual ligand. Further 
impurities were removed from the filtered product by Soxhlet 
extraction (ethanol). The purified product was then dried at 
150 ° for 24 h in vacuo. Less defective samples were prepared 
using larger concentrations of terephthalic acid: 3.0, 4.0, 4.5 
and 5.0 mmol. All samples were loaded in 0.5 mm quartz glass 
capillaries with AI 2 O 3 as an internal temperature standard (ca 
50% by mass) and sealed with cotton wool. 

Thermogravimetric Analysis 

Thermogravimetric analysis was carried using a Mettler 
Toledo TGA/DSC 1 System, heating from 30 °C to 650 °C at a 
rate of 20 °C min -1 before cooling down to room temperature. 

Variable-temperature synchrotron X-ray diffraction mea¬ 
surements 

Variable temperature X-ray powder diffraction measurements 
were carried out on the Ill beamline at the Diamond Light 
Source (k = 0.82562 A) using a position sensitive detector 
with a Cyberstar hot air blower .^ 1 52 Analysis of the data was 
carried out using Topas Academic version 4.1 . 53 All samples’ 
purity were checked using a Rietveld fit to the reported struc¬ 
ture for UiO-66(Zr)E3 w ith Hf substituted for Zr. The addi¬ 
tional primitive peaks were fitted using a secondary Pawley 
phase with the same lattice parameter but with additional size 
broadening as described in Ref. 6 Lattice parameters were 
determined using Pawley refinement in space group Fm3m , 
omitting the low angle region where the broad primitive reflec¬ 
tions make an important contribution. The Pawley refinements 
were carried out using the previous temperature’s refined pa¬ 
rameters as the input for the next refinement. 

Quantum Mechanical Calculations 

All structures were fully relaxed by optimizing both atomic 
positions and unit cell parameters. We performed quantum 
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mechanical calculations in the density functional theory ap¬ 
proach with localized basis sets (CRYSTAL14 code 55 ). We 
used the B3LYP hybrid exchange-correlation functional 56 and 
all-electron basis sets for all atoms,^ 5 except Hf for which we 
used a small relativistic effective core potential with 12 outer 
electrons considered explicitly. 57 Full methodological details 
can be found in in Refs. [6] and I57J The optimised structures 
are included as crystallographic information files as part of the 
ESIf. 
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